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Abstract 

We present a systematic analysis of the three-flavor Mikheyev-Smirnov- 
Wolfenstein (MSW) oscillation solutions to the solar neutrino problem, in 
the hypothesis that the two independent neutrino square mass differences, 
5m 2 and m 2 , are well separated: 5m 2 <C m 2 . At zeroth order in 5m 2 /m 2 , 
the relevant variables for solar neutrinos are 5m 2 and two mixing angles, oj 
and <f>. We introduce new graphical representations of the parameter space 
(5m 2 , u, 4>), that prove useful both to analyze the properties of the electron- 
neutrino survival probability and to present the results of the analysis of 
solar neutrino data. We make a detailed comparison between the theoretical 
predictions of the Bahcall-Pinsonneault standard solar model and the cur- 
rent experimental results on solar neutrino rates, and discuss thoroughly the 
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MSW solutions found by spanning the whole three-flavor space (5m 2 , u>, <j>). 
The allowed regions can be radically different from the usual "small mixing" 
and "large mixing" solutions, characteristic of the usual two-generation MSW 
approach. We also discuss the link between these results and the indepen- 
dent information on neutrino masses and mixings coming from accelerator 

and reactor oscillation searches. 

PACS number(s): 26.65.+t,14.60.Pq,13.15.+g. 



2 



Typeset using REVTf^X 



I. INTRODUCTION 



The measurement of the Z — > vv width at the CERN Large Electron- Positron collider 
[ has shown conclusively that there are three generations of light neutrinos. However, 
the problem as to whether the neutrinos have non-zero mass and mixing, as is the case for 
quarks, remains one of the most prominent experimental and theoretical problems in particle 
physics. 

On the one hand, there is a vast and variegated set of experimental information, ranging 
from the negative results of laboratory neutrino oscillation searches at accelerators and 
reactors, to the possible indications for flavor transition processes in natural beams such as 
solar and atmospheric neutrinos. On the other hand, current theoretical models (e.g., the 
see-saw mechanism ||) allow many different textures for the three-generation neutrino mass 
matrix. Thus, it seems appropriate, at this stage, to keep to a minimum the number of a 
priori assumptions in the interpretation of the experimental results. 

In this work, we focus on the three-flavor Mikheyev-Smirnov-Wolfenstein (MSW) || 
neutrino oscillation solutions to the solar neutrino problem [|J. The only hypothesis that we 
make is to assume a neutrino mass spectrum as in Fig. 1. More precisely, we assume one of 
the two independent square mass differences, m 2 , to be very large (more than one order of 
magnitude higher) as compared to the smallest one, 5m 2 , which drives the matter-enhanced 
oscillations. 

At zeroth order in 5m 2 /m 2 , the relevant MSW parameters (see Sec. Ill) are 5m 2 and 
two mixing angles, u and 0, ranging in the first quadrant [0, vr/2]. We introduce new 
graphical representations of the parameter space (5m 2 , uj, <fi), that prove useful in analyzing 
the properties of the three-flavor MSW probability, as well as in showing the detailed results 
of fits to the experimental data. The parameter space is studied exhaustively, and all the 
solutions are discussed thoroughly. We find three-flavor solutions to the solar neutrino 
problem that are considerably different from the usual two-flavor results. 

The present study is more complete and detailed than previous three-flavor phenomeno- 
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logical analyses of the solar neutrino problem performed by other authors |5|-|9| and by 
ourselves [p)|-|T2|j. The work presented in this paper is part of a wider research program, 
in which we intend to analyze the world neutrino oscillation data under the sole hypothesis 
represented by the spectrum in Fig. 1. The analysis of the most constraining accelerator 
and reactor experiments has already been completed ]nj. We discuss the link between the 
information obtained in the present work from solar neutrinos and the results obtained in 
T3fl from laboratory oscillation searches. Concerning atmospheric neutrinos, a complete 



three-flavor analysis (partially addressed in |iOHl4f) is in progress. 

This paper has the following structure. In Sec. II we review the current theoretical expec- 
tations and experimental observations of the solar neutrino fluxes, and show for completeness 
the usual two-generation MSW fits to the data. In Sec. Ill we present the three-flavor oscil- 
lation formalism and show the new graphical representations of the parameter space. The 
properties of the v e survival probability are also discussed. In Sec. IV we explore the param- 
eter space exhaustively, and find the three-flavor solutions to the solar neutrino problem. 
Both in Sec. Ill and in Sec. IV we show how solar neutrino results can be compared to the 
information coming from accelerator and reactor neutrino experiments. We conclude and 
summarize our work in Sec. V. 



II. SOLAR NEUTRINO PROBLEM AND TWO-GENERATION MSW 

SOLUTIONS 

In this section we present the experimental data used in this work, and compare them 
with the theoretical expectations. Particular attention is paid to correlation effects. Al- 
though our paper is focused on the three-flavor MSW mechanism, we also show, for the 
sake of completeness, our fit to the data in the simpler and more familiar two-generation 
approach. 
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A. Experimental data and theoretical predictions 



There are four operating solar neutrino experiments: The Homestake [15] chlorine (CI) 



detector, the GALLEX [16] and SAGE 17] gallium (Ga) detectors, and the Kamiokande |H 
water- Cherenkov detector. The observed solar neutrino rates are shown in Table I, along 
with the corresponding theoretical values as predicted in the recent refined Standard Solar 
Model (SSM) by Bahcall and Pinsonneault (BP) [|HJ. The significant deficit of observed 
neutrinos constitutes the well-known "solar neutrino problem" M. 

The information contained in Table I does not characterize completely the significance 
of the solar neutrino problem, since the theoretical uncertainties (last column of Table I) 
are not independent. It has been shown in a previous work that the correlation of the 
uncertainties can be calculated analytically, starting from the errors affecting the SSM input 
parameters and the neutrino capture cross-sections. 

In Fig. 2 we show the experimental and theoretical 99% C.L. contours (A% 2 = 9.21, 
Ndf = 2) for any couple of experiments (the GALLEX and SAGE data have been combined 
quadratically) . The projections of the ellipses onto the coordinate axes represent, with good 
approximation, the ±3<r errors for each experiment separately. Notice the strong correlation 
of the theoretical errors. If the correlations were neglected, the theoretically allowed regions 



would be larger and the discrepancy with the experiments would be underestimated 



The last relevant piece of data is the night-day (N — D) asymmetry of the solar neutrino 
flux as measured in the Kamiokande detector ||23JT8|1: 



N - D 

N + D = 0.07 ± 0.07(stat) ± 0.04(syst) . (1) 

The above value is consistent with no asymmetry (i.e., no oscillations), and thus contributes 
to exclude a part of the neutrino oscillation parameter space. 
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B. Two-generation MSW analysis 



The MSW mechanism [EJ of matter-enhanced neutrino oscillations is one of the most 



promising candidates for the solution of the solar neutrino problem (for reviews, see [25 



In the hypothesis that only two neutrino generations are involved in the v e transitions, the 
relevant parameters are one mass square difference, 5m 2 , and one mixing angle, uj. 

Even in the simple two-family scenario, the calculation of the electron neutrino survival 
probability, P{y e —>■ u e ), as well as the data analysis, are not trivial tasks and require refined 
computer codes. We take this opportunity to mention a few state-of-the-art features of our 
codes. We smear P{y e — > u e ) over the neutrino production region in the most general way, 
i.e., we also consider neutrino production points off the line that joins the center of the sun 
and the detector. This requires a double (radial and azimuthal) space integration in the sun, 
as well as the evaluation of the solar electron density gradient along off-radial trajectories 
(instead of the simple derivative). The possibility of a double resonance (see, e.g., |2^j2^)) is 
included. We also include the (computer time-consuming) earth regeneration effect at 



each detector location, and in particular we calculate the night-day asymmetry |3T[ a t the 
Kamiokande site. The earth electron density is modeled as a step function with five steps 
|I(J , corresponding to the relevant radial shells [5^] . The theoretical uncertainties and their 
correlation effects are taken into account in the statistical x 2 -analysis of the data | |20| | . 

The results of our two-flavor analysis are shown in Fig. 3. The thick, solid lines represent 
the contours of the regions allowed at 95% C.L. (A% 2 = 5.99, Ndf = 2). In the first 
four panels, dotted lines represent curves of iso-signal. In the last two panels, all data are 
combined without and with theoretical errors, and the allowed regions at 90%, 95% and 
99% C.L. are shown. Notice the widening of the "small" and "large" mixing angle solutions 
when the theoretical uncertainties are taken into account. A third solution appears at 99% 
C.L. for large angles and small 8m 2 . 

When all the uncertainties are included, the y 2 function reaches its minimum, Xmm = 
0.74, at sin 2 2a; = 8.06 x 10~ 3 and 5m 2 = 5.21 x 10~ 6 eV 2 . The secondary (large angle) 
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minimum is reached at sin 2 2u = 0.64 and 5m 2 = 1.45 x 10~ 5 eV 2 , with a value Xmm = 1-89— 
not a bad fit. 

The MSW solutions in Fig. 3 compare well with the most recent two-generation analyses 



of the solar neutrino problem p2| , p3| , p3| -|38[] , modulo small differences due to the input data 
and their treatment. In particular, we have verified that, excluding the earth regeneration 
effect, our results become very similar to those obtained very recently by Bahcall and Krastev 
in [381, where the same solar model and experimental data were used as input. 



III. THREE-FLAVOR MSW SOLAR NEUTRINO OSCILLATIONS: 

THEORETICAL ASPECTS 

In this Section we recall the formalism that will be used in the three-flavor analysis, 
under the hypothesis that the neutrino mass spectrum is as in Fig. 1. We introduce new 
graphical representations of the parameter space, and use them to discuss the properties of 
the three-flavor MSW probability P[y e — > u e ). 



A. A simple three-flavor framework and its parameter space 



In the general three-generation case, the flavor eigenstates u a (a = e, fi, r) are a super- 
position of three mass eigenstates z/j {i = 1, 2, 3): v a = \J ai Vi. The (unitary) neutrino mixing 
matrix U a i is usually parameterized in the same way as the standard Cabibbo-Kobayashi- 



Maskawa (CKM) mixing matrix in the quark sector | 39|[ , involving three mixing angles, 
(#12, #23, #13) G [0, tt/2], and one CP-violating phase, 5. Redefining the angles as uj = 6*12, 
•0 = #23, and 4> = #13, the standard parameterization reads: 



( 



\ 



COS UJ COS ( 



sin uj cos ( 



— sin uj cos if) — cos uj sin tp sin (f>e lS cos uj cos tp — sin uj sin ip sin 4>e lS sin cos 4> 
sin uj sin ip — cos uj cos if) sin <f>e lS — cos uj sin if) — sin uj cos if) sin <f>e zS cos ^1 cos (f) 



( \ 



(2) 
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Permutations of the mass eigenstate labels (1, 2, 3) do not change the physics, provided 
that the mixing angles (uj, ip, <fi) are left to range in the first quadrant |[39|| . The parameter 
set is completed by two independent neutrino square mass differences, that we write as: 

5m 2 = m\ — m\ , m 2 = m 2 — m\ , (3) 

so that the general three-flavor parameter space is (5m 2 , m 2 , uj, ip, <fi, 5). 

In principle, a complete phenomenological analysis of three-flavor neutrino oscillations 
should span the entire 6-dimensional manifold (5m 2 , m 2 , uj, ip, <fi, 5). This task would be 
exceedingly difficult — and perhaps not really useful — at this time, due to vastness of the 
parameter space, as opposed to the scarcity of evidences in favor of neutrino oscillations. 

A much more manageable framework is obtained under the simple hypothesis that the 
mass spectrum has the property shown in Fig. 1: 5m 2 <C m 2 . The two cases (a) and 
(b) in Fig. 1 lead to the same solar neutrino physics. We will always refer to (a) in the 
following, and conventionally label the mass eigenstates in the order of increasing mass 
(mi < m 2 < m 3 ). In this far as solar neutrinos are concerned, the two quasi- 

degenerate mass eigenstates v\ and v 2 participate actively to matter-enhanced transitions, 
while the "lone" state v 3 is a spectator neutrino that reveals its presence only through the 
mixing |^0|-^2] . 



At zeroth order in 5m 2 jm 2 , the following properties hold (see [|T(],[ni and references 
therein): (1) solar neutrino experiments probe only the subspace (5m 2 , uj, 0); (2) terres- 
trial (i.e., accelerator, reactor, atmospheric) neutrino experiments probe only the subspace 
(m 2 , ip, 0); (3) effects related to the CP-violating phase 5 are unobservable, so that the 
mixing matrix in Eq. (2) can be taken as real for our purposes. 

In other words, solar neutrino experiments explore the space (5m 2 , U 2 X , U 2 2 , U 2 3 ) that 
is the small mass square difference and the three matrix elements related to v e (the fast 
m 2 -driven oscillations being effectively averaged out) ||10| . Terrestrial (accelerator, reactor, 
atmospheric) neutrino experiments explore instead the space (m 2 , U 2 3 , U 2 3 , U 2 3 ) that is the 
large mass square difference and the three matrix elements related to v 3 (the slow 5m 2 -driven 



oscillations being effectively frozen) [pT3| . 

The five mixing matrix elements U e i and probed by solar and terrestrial oscillation 
searches satisfy the two unitarity conditions: 

U^ + U% + U^ = l , (4a) 

U 2 el + U 2 e2 + U 2 e3 = l . (4b) 

In [[b| it has been shown that Eq. (4a) can be embedded in a triangular representation of 
the space spanned by the the "lone" neutrino state: u 3 = U e3 Ve+U^ i^^+U t3 v t . In Fig. 4, we 
embed analogously the second unitarity condition, Eq. (4b), in a triangular representation of 
the space spanned by the electron neutrino state: v e = U e \ v\ + U e2 v 2 + U e3 u 3 . The triangle 
in Fig. 4 has equal sides, unit height, and corners corresponding to the mass eigenstates v%, 
z/ 2 , z/ 3 . A generic state v e is represented by a point in the triangle. The are identified with 
the heights projected from v e . The sum of such heights is always equal to the total (unit) 
height of the triangle. In Fig. 4 we also chart, in the lower subplot, the triangle coordinates 
in terms of uo and 0. The usual two-generation limit is reached for = 0, that is for v e on 
the side joining v\ to v% (and thus decoupled from z/ 3 ). 

In Fig. 5 we show a synoptic presentation of the two triangular graphs introduced in Fig. 4 
of this work and in Fig. 2 of fi~3fl . On the left (right) side we display the parameterization of 
solar (terrestrial) neutrino oscillations, including the expansion of the relevant mixing matrix 
elements in terms of the mixing angles. Note that both solar and terrestrial neutrinos probe 



the element U e s, i.e., the angle 0, as discussed in [flO| p~3|1 and also in |43| , |44]| . This important 
link between the experiments sensitive to 5m 2 (solar) and those sensitive to m 2 (terrestrial) 
will be utilized in Sec. IV. 



B. Electron neutrino survival probability 



The electron neutrino survival probability, P[y e — > z/ e ), is the fundamental quantity to 
compute in the solar neutrino analysis. In this section we study some properties of the three- 
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flavor MSW survival probability, P^, sw , as a function of u, 0, and Sm 2 /E u (E„ being the 
neutrino energy). We start with a brief review of a well-known analytical approximation to 
Pj^ sw , that has been used in the present work. In order to avoid unnecessary complications, 
in this section and in the related figures we ignore temporarily the smearing of -Pj^ sw over 
the neutrino production region, and the earth regeneration effect. These effects will be 
included, however, in the global analysis of Sec. IV. 

In the limit 5m 2 <C m 2 , the probability -P^ sw takes the simple form []40,5|: 



pMSW = cos 4 pMSW + sin 4 5 (5) 

where P^, sw is the MSW probability in the two-generation limit, provided that the solar 
electron density N e (x) at any point x is effectively replaced by N e (x) cos 2 0. 

A very accurate and widely used analytical approximation to P^, sw * s the so-called 



"Parke's formula" 45 



f^ sw = i+fi-eP c )cos2wcos2o;S l . (6) 



2 V2 

This approximation was also studied by other authors [4"B|. In Eq. (6), Pq is the Landau 



Zener-Stueckelberg crossing probability |^7|, that we use in the improved form f4"5| , f4T| valid 
for a close-to-exponential radial density; oj^ is the mixing angle in matter at the production 
point. The step function [^J switches from (non- resonant propagation case) to 1 
(resonant propagation case) if the mixing angle uj in matter can assume the value 7r/4 at a 
point x along the neutrino trajectory in the sun: 

w m (x) = 7 « 5m 2 ^^- = 2V2G F N e (x)E„ . (7) 
4 cos^ (p 

Equation (7) numerically reads 

5m 2 cos2cj - E u NJx) , , 

where N e (0) is the electron density at the center of the sun [ |T9|j . 

A technical remark is in order. In the literature, the numerical validity of the approxi- 
mations implicit in Eq. (6) has been discussed |49|j50[1 mainly for uj (and <p P7| ) in the first 
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octant [0, 7r/4]. Since we are interested in analyzing the full mixing angle space (the first 
quadrant), we have calculated the function Pj)J = P 3 ^, sw (5m 2 / E^) with a Runge-Kutta 
integration of the three-flavor MSW equations, for a number of representative (u;, 0) val- 
ues in [0, 7r/2]. The agreement between the analytical approximation of Eqs. (5,6) and the 
"exact" (but much more time-consuming) Runge-Kutta calculation is as good in the second 
octant as it is in the first, provided that 5m 2 / ' E u is taken above ~ 10 -8 eV 2 /MeV . For 
values of 5m 2 /E v lower than ~ 1CT 8 eV 2 /MeV, one should more properly consider just-so 



vacuum oscillations 51 



There are a few interesting limits for P 3 ^ sw . If 5m 2 jE v > 1(T 5 eV 2 /MeV, the MSW 
mechanism is not effective and P^r tends to its vacuum value P 3 ™ c given by 

p™ c = cos 4 P™ c + sin 4 

= cos 4 <j) (l - ^ sin 2 2uj + sin 4 (9) 
= 1 - 2(U 2 el U 2 2 + U 2 2 U 2 3 + U 2 3 U 2 el ) , 

where the fast (5m 2 / E^-driven oscillations have been averaged to 1/2. 

Given the above expression for P 3 ™ c , and Eqs. (5,6) for P 3 ^ sw , it follows that P|/! sw = 
P™, for arbitrary values of (5m 2 /E u ), in (at least) three cases: 

u = — ► P 3 ^ SW = cos 4 + sin 4 = P™ C , (10a) 
*=\ — P 3 ^ SW = ^cos 4 + sin 4 = P 3 r , (10b) 
oo = ^ — > P 3 f W = cos 4 + sin 4 = P™ C . (10c) 

The above properties are exact, i.e., they do not depend on the approximations implicit 
in Eq. (6). A subcase of Eq. (10b) is obtained in the "maximal mixing scenario" ||52|| . 
corresponding to the center of the triangle in Fig. 4 (sin 2 uj = 1/2, sin 2 0= 1/3). 
A further property follows from Eqs. (5,6,9): 

1 



MSW 



cos 4 + sin 4 = P rac 



J rn 4 4 

This property, however, is not exact and has the same limits of validity of Eq. (6) (see also 
the discussion at the end of this section). 
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The equations (9-11) are helpful in understanding the behavior of P^ sw as a function 
of its arguments (Sm 2 /E u , uj, <fi), as we discuss now in some detail. 

In Fig. 6, curves of constant Pj^ sw are shown for selected values of 5m 2 /E u in the 
triangular representation. In the first subplot (Sm 2 /E u ^> 10~ 5 eV 2 /MeV), it is -Pj^ sw = 
P 3 ™° in the whole triangle, and the iso-probability curves are circles [as it is easily derived 
from Eq. (9)], with a minimum value P^ u = 1/3 at the triangle center. In the second 
subplot (5m 2 /E u = 1.75 x 10 -5 eV 2 /MeV), the curves begin to be deformed by the MSW 
effect for v e close to v\ (the low left corner). Notice, however, that just along the lines 
at uj = 0, 7r/4, 7r/2 (refer to Fig. 4 also) the probability takes the same values as in the 
previous plot (-P 3 ^ sw = P™), as expected from Eqs. (lOa-c). These properties also hold 
for the remaining four subplots where, however, the increasingly important MSW effect 
prevents a graphical resolution of the iso-lines at small uj values, i.e., close to the left side of 
the triangle. 

In Fig. 7 we map the triangular parameter space onto a "square" bilogarithmic plot with 
coordinates (tan 2 u, tan 2 0). This new representation has the advantage that the "triangle 
corners" are infinitely expanded, although the nice symmetry properties of the triangular 
representation are lost. A similar representation was introduced in |T\| to map the parameter 



space of terrestrial neutrinos. The values of 5m 2 / E v and the iso-probability curves in Fig. 7 
correspond exactly to those in Fig. 6. In Fig. 7 we also draw two curves (thick, dashed lines) 
corresponding to uj^ = 7r/4 [Eq. (7) with x = 0] and uj = ir/4. These curves separate the 
zones where the propagation is non- resonant from those where it is resonant, as indicated 
in all panels. The behavior of P^, is most interesting near the curve at uj^ = ir/4 and in 
the resonance region. Notice, for the subplots at 5m 2 / E v = 1.45 x 10~ 5 eV 2 /MeV or lower, 
the appearance of a depletion zone for P|^ sw at small values of uj (the "small angle MSW 
solution"). 

In Fig. 7, the properties expressed by Eqs. (lOa-b) can be easily checked. The implications 
of Eq. (11) are instead more subtle; if one takes any two points at the same height (same 
0) on the dashed curves {uj = ir/4 and uj^ = vr/4), then the probabilities P^ v at these 
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conjugate points are equal. It may be said, in a figurative way, that the MSW mechanism 
"cuts" the (u, 4>) plane along the dashed lines at to = 7r/4 and oJ^ = 7r/4, "separates" 
the edges (that would coincide for vacuum oscillations) and "fills" the cut zone with a new 
resonant region where the survival probability can be much lower than in the vacuum case. 
The probabilities in the "old" regions at the right and left of the two dashed lines are only 
weakly (adiabatically) different from the corresponding vacuum oscillation values. 

We finally show in Fig. 8 the behavior of -P 3 ^ sw for six representative values of 0, in a 
bilogarithmic (Sm 2 /E 1/ , tan 2 u) mass-mixing plane (a perhaps more familiar representation). 
As in the previous figure, the two thick, dashed lines correspond to the curves at uj = tt/4 
and uP m = 7r/4, that separate non- resonant and resonant propagation regions. The curves of 
iso-probability in the resonant region have the characteristic triangular shape f26|. Notice 



that the variations of -Pj^ sw within the mass-mixing plane decrease as increases and the 
MSW effect is suppressed. 



A final remark about the calculation of the survival probability. As shown in [^9 
Eq. (6) is a good approximation except for the particular case of a neutrino created close 
to the resonance (u^ ~ tt/4) and with very small vacuum mixing (tan 2 cu < 10~ 4 ). This 
situation would correspond to a very narrow strip, localized along the thick, dashed curves 
in Figs. 7 and 8 for tan 2 c<j < 10~ 4 . The local discrepancy with the exact probability is, 
however, unimportant in the calculation of the neutrino rates, being effectively suppressed 
by the integration over the neutrino production region and energy spectrum, and being 
confined to values of u where no solution to the solar neutrino problem is found. 

IV. THREE-FLAVOR MSW SOLAR NEUTRINO OSCILLATIONS: 

PHENOMENOLOGY 



In this section we present the results of our global three-flavor MSW analysis of solar 
neutrino data (Table I) in the parameter space (5m 2 , uj, <fi). We span the coordinates 5m 2 ® 
tan 2 uo ® tan 2 in the range [10 -7 eV 2 , 10" 3 eV 2 ] <g> [1(T 6 , 1(T 2 ] ® [1(T 4 , 10~ 4 ] and analyze 
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the data through a x 2 statistic, including experimental and theoretical errors and their 
correlations. 

We find the absolute minimum, Xm in = °- 74 > at $ m2 = 5 -21 x 10~ 6 eV 2 , u = 2.57°, = 0°. 
It coincides with the minimum of the two-generation analysis (Sec. II B). The surfaces at 
x 2 = x 2 m . m + Ax 2 , with Ax 2 = 6.25, 7.82, and 11.34, define the regions allowed at 90%, 
95% and 99% CL. for three degrees of freedom. We present representative sections of these 
three-dimensional allowed regions in the bilogarithmic mass-mixing plane (5m 2 , tan 2 co>) at 
several values of tan 2 0, and in the bilogarithmic mixing-mixing plane (tan 2 u, tan 2 0) at 
several values of 5m 2 . We also show representative fits to separate pieces of data. We 
finally discuss how this solar neutrino analysis is linked to the terrestrial neutrino oscillation 



constraints obtained in 13 



A. Analysis in the mass-mixing plane 

In Fig. 9 we show the sections of the allowed regions in the bilogarithmic plane 
(tan 2 a;, 5m 2 ), for representative values of tan 2 0. Notice than both u and are allowed 
to range above the first octant in the fit. In the first panel, tan 2 is zero and purely two- 
generation oscillations take place. In fact, this subplot is a mapping of Fig. 3, modulo the 
different number of degrees of freedom in the definition of the CL. contours. A moderate 
increase of tan 2 from 0.1 up to 0.4 produces an increasing deformation of the allowed 
contours. The contours tend to merge in the upper part, corresponding to the upper side 
of the MSW "triangles" of Fig. 8. The two separate solutions merge in a single connected 
solution at tan 2 ~ 0.5, as also shown in |7|,(|. For increasing value of tan 2 0, the allowed 
region becomes broader and less structured, and then rapidly shrinks and disappears at 
tan 2 > 1.5. Figure 8 helps in understanding the flattening of the allowed region in Fig. 9 
at the higher values of tan 2 0. 

From Fig. 9 it can be noticed that, although there are no solutions for u > 7r/4, the 
fit can still be acceptable for > 7r/4. In Fig. 10 we then show the values of x 2 ~ Xmin 
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as a function of tan 2 0, for 5m? and tan 2 uj unconstrained. The minimum is reached at 
= 0°, meaning that the two-generation limit is preferred by the present data. However, 
three-flavor solutions are allowed for values of tan 2 as high as 1.4 (i.e., — 50°) at 90% 
CL. There are no acceptable solutions at 99% CL. for tan 2 0^2. 

The behavior of the solutions in Fig. 9 can be understood better by separating the 
information coming from single pieces of data. In Fig. 11 we show the separate fits to the 
gallium, chlorine and Kamiokande (total rate and night-day asymmetry) measurements, in 
the same coordinates as in Fig. 9, for three values of tan 2 0. The thick, solid line represents 
contours at 95% CL. The allowed regions are marked by stars. The dotted lines represent 
contours of iso-signal. The structured shape of the iso-signal curves is globally similar to the 
iso-probability curves in Fig. 8. The regions allowed by the separate pieces of data become 
wider for increasing 0, and even extend above u = tt/4, but their mutual compatibility do 
not necessarily increase. In particular, notice in Fig. 11 that the chlorine data play a major 
role in confining u below tt/4 in the global fit of Fig. 9. 

We have found three-flavor solutions to the solar neutrino problem at relatively large 0. 
These solutions are characterized by a survival probability, Pj^ sw , that varies more slowly 
with energy as compared to the two-generation case (see Fig. 8). Probing the large-0 three- 
flavor solutions will therefore represent a formidable challenge to the new generation of solar 
neutrino experiments, SuperKamiokande f53 |, the Sudbury Neutrino Observatory (SNO) 
54| , and the Imaging of Cosmic and Rare Underground Signals (ICARUS) f55| , that are 



meant to probe the energy-dependence of the solar neutrino flux at earth. 



B. Analysis in the mixing-mixing plane 

In Fig. 12 we show the sections of the allowed region in the mixing-mixing plane 
(tan 2 a;, tan 2 0) for representative values of 8m 2 . In all panels, the two-generation limit is 
reached for tan 2 0^0 (lower side), and the solutions lie along one or more of a three-sided 
ri-shaped band. The sides correspond to "small u" "large u" and "large 0" solutions. The 
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"small 00" and "large oj n solutions protrude in the three-flavor space from the usual MSW 
solutions in the two-generation limit. The "large (ft" solution is genuinely three-flavor and 
corresponds to the horizontal depletion zone of the survival probability Pj^ sw in Fig. 7. In 
the first two and last two subplots there is no solution in the two-generation limit. 

In Fig. 13 we show the fits to single experimental data, in the same coordinates as in 
Fig. 12, for three values of 8m 2 . As in Fig. 11, the thick, solid lines correspond to iso-Ax 2 
contours at 95% C.L., and the dotted lines represent iso-signal curves. The allowed regions 
are marked by stars. Notice how, in all cases, the solutions in the two-generation limit (lower 
side, (j) — > 0) merge as (ft increases. 



C. Solar neutrinos vs accelerator and reactor neutrinos 

In this work we have found and discussed the MSW solutions to the solar neutrino 
problem in the parameter space (5m 2 , tan 2 u;, tan 2 0). In the work jHJ, the results of the 
established accelerator and reactor neutrino oscillation searches were analyzed in order to 
constrain the parameter space (m 2 , tan 2 ^, tan 2 0). Both works are based on one, and only 
one, hypothesis about the neutrino spectrum, namely 5m 2 <C m 2 (as in Fig. 1). The angle (ft 
is common to both spaces (8m 2 , tan 2 tu, tan 2 (ft) and (m 2 , tan 2 ?/>, tan 2 0), as also evidenced 
in Fig. 5. 



in 



We invite the reader to place, side-by-side, Fig. 12 of this work and Fig. 10 of |13 
the way suggested by Fig. 5. These two figures provide us with 12 x 12 combinations of 
(8m 2 , m 2 ). For any given couple of square mass differences (8m 2 , m 2 ), one can find in these 
two figures the allowed regions in the (tan 2 a;, tan 2 (ft) and (tan 2 ?/>, tan 2 (ft) mixing-mixing 
planes. We have purposely chosen the same vertical scale (same decades in log tan 2 (ft) in 
both cases. The obvious condition that tan 2 (ft (the vertical coordinate) must be the same in 
both figures imposes a strong compatibility constraint. This constraint eliminates some of 
the solutions that are individually allowed by solar or terrestrial neutrino oscillation searches. 
For instance, the regions allowed by accelerator and reactor data at very large tan 2 (ft (i.e., 
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— > 7r/2) are never allowed by solar neutrino data, that keep tan 2 below ~ 1.5. In other 
words, if m 2 is in the range of sensitivity of accelerator and reactor oscillation searches, then 
the comparison of terrestrial and solar data requires relatively small values of 0. Another 
example is that the MSW-allowed values Sm 2 ~ 1.5 x 10~ 4 eV 2 and 8m 2 ~ (2.2-3.2) x 1CT 6 
eV 2 are not compatible with accelerator and reactor data for m 2 > 10~ 2 eV 2 . One can 
explore various other (in) compatibility situations for the 12 x 12 cases of (8m 2 , m 2 ) that 
can be extracted from the figures discussed above. 

We recall that the only hypothesis adopted in our analysis is Sm 2 <C m 2 . Virtually, all 
see-saw models of neutrino mass generation satisfy this condition. Model builders can thus 
check if their favorite values of (Sm 2 , m 2 , uo, ip, 0) fall in the allowed regions determined 



in this work and in [13|. Here we discuss briefly the "threefold maximal mixing" model, 



recently investigated in [56 



The maximal mixing scenario is realized in the centers of the triangles in Fig. 5, that 
is for [T 2 = U 2 3 = 1/3 or, equivalently, for (tan 2 u;, tan 2 ?/>, tan 2 0) = (1, 1, 1/2). Figure 9 
shows that the combination (tan 2 uo, tan 2 0) = (1, 1/2) is excluded at more than 99% C.L. by 
solar v data (within the MSW approach). The maximal mixing scenario could be reconciled 
with solar v data only by dropping the Homestake measurement, that keeps u below 7r/4. 
From Fig. 10 in ||13|| , it follows that accelerator and reactor data exclude the maximal mixing 
combination (tan 2 ip, tan 2 0) = (1, 1/2) for all values of m 2 in their present sensitivity range 
(m 2 > 10~ 2 eV 2 ). Thus, the maximal mixing is not supported by all solar neutrino data, and 
is (trivially) compatible with accelerator and reactor experiments only below their present 
m 2 -sensitivity limits. 



V. SUMMARY AND CONCLUSIONS 

We have discussed the three-flavor MSW mechanism, in the hypothesis 6m 2 -C m 2 and for 
mixing angles uo and in their maximal range [0, vr/2]. New useful graphical representations 
of the parameter space (Figs. 4-6) have been introduced. A complete analysis of present 
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data in the mass-mixing space (6m 2 , uo, 0) has been performed, and the solutions have been 
thoroughly discussed. The data favor the two-generation limit (0 = 0°), but a wide range 
of three-flavor solutions are allowed up to — 50°. The main results are summarized in 
Figs. 9 and 12. 

The u e oscillation probability has been studied as a function of the mass and mixing 
parameters (Figs. 7 and 8). In the allowed regions at large 0, the three-flavor probability 
varies more slowly with energy than in the two-generation case. Thus, it will be difficult to 
study the large-0 solutions in the future solar neutrino experiments that are meant to probe 
the energy dependence of the solar neutrino flux. 

It has been shown that the solar neutrino results are linked to those obtained in an 
analysis of accelerator and reactor neutrino oscillation searches |13[], where the mass-mixing 
space (m 2 , ip, 0) has been thoroughly explored (Fig. 5). This link can be used to further 
constrain the allowed ranges of the parameters. In particular, for m 2 > 10~ 2 eV 2 , the 
angle is constrained to relatively small values. The threefold maximal mixing scenario 
pEfl is disfavored by solar neutrino data (in the whole explored range of 8m 2 ) as well as by 
accelerator and reactor experiments for m 2 > 10~ 2 . 

The present study and the work in [TB[] represent a detailed and accurate analysis of 
solar, accelerator, and reactor neutrino oscillations in three-flavors. The analysis has made it 
possible to obtain significant bounds in the mass-mixing parameter space (5m 2 , m 2 , u, ip, 0). 

When this work was completed, we become aware of a recent three-flavor analysis of 



solar neutrino data performed by Narayan et al. in ||57|| ■ We agree only qualitatively on 
their results about the MSW solutions. We do not agree on the precise shape of the allowed 
regions, and do not confirm their claim that the solutions are dramatically enlarged when 
the confidence level is moderately increased (from la to 1.6a in their work). The information 
contained in []571 is not detailed enough to trace the differences. 
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TABLES 



TABLE I. Neutrino rates observed in the four operating solar neutrino experiments [15 18], 
and corresponding theoretical predictions of the Standard Solar Model (SSM) by Bahcall and 
Pinsonneault [Oil. 



Experiment 



Ref. 



Measured rate ±(stat)±(syst) 



Units 



SSM rate @ 



Homestake 

GALLEX 

SAGE 

Kamiokande 



m 



2.55 ±0.17 ±0.18 



77.1 ±8.51^4 



69 ± llil 



2.89 ±0.22 ±0.35 



SNU 
SNU 
SNU 
10 6 cm-V 1 



9.3 
137 



1.2 
1.4 

+8 



1371? 



6.62; 



-0.93 
-1.13 
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v mass spectrum hypothesis! 6m 2 « m 2 
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FIG. 1. The neutrino mass spectrum assumed in this work. The hypothesis 5m 2 <C m 2 is 
independent of the zero of the absolute mass scale. Solar neutrinos do not distinguish the two 
cases (a) and (b) at zeroth order in 5m? /m 2 . 
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Solar neutrino problem(s) 




99 % C.L. contours 
(2 DF in each plane) 
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FIG. 2. The regions allowed at 99% C.L. by the present solar neutrino data (dashed lines) and 
by the standard solar model (SSM) of Bahcall-Pinsonneault |TI| (solid lines). The coordinates are 
the chlorine (CI), gallium (Ga), and water-Cherenkov (normalized Kam/SSM) signals. 
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FIG. 3. Two-generation MSW solutions to the solar neutrino problem. In the upper four panels 
(fits to single data), the solid lines define the regions allowed at 95% C.L. for two degrees of freedom 
(including theoretical errors), and the dotted lines represent iso-signal contours. In the lower two 
panels (fits to all data), the MSW solutions are shown at 90%, 95%, and 99% C.L., without and 
with theoretical errors. 
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FIG. 4. Triangular plot representing the neutrino state v e in terms of its massive components 
v \i v 2i v-&- For a triangle of unit height, the unitarity relation + C/ 2 2 + U% = 1 is enforced 
by identifying the with the three distances of v e from the sides (dashed lines). The usual 
two-generation oscillation limit is reached for v e onto the lower side joining v\ to vi- Iso- lines of 
sin 2 a; and sin 2 <fi (parameterizing the elements C/ 2 j) are charted in the lower plot. 
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FIG. 5. A synoptic presentation of the parameter space for solar and terrestrial (accelerator, 
reactor, atmospheric) neutrino oscillation searches. Solar neutrinos probe the small square mass 
difference, 5m 2 , and the mixing matrix elements U^. Terrestrial neutrinos probe the large square 
mass difference, to 2 , and the mixing matrix elements U 2 3 . Notice that f/ 2 3 is probed in both cases, 
whilst the CP-violating phase 5 is unobservable. The expansion of in terms of the mixing 



angles (uj, <j), ifj) is shown. The reader is referred to the text and to 13] for further details. 
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3v MSW = P{v t ^v t ) 





FIG. 6. Iso-lines of survival probability P{y e — > f e ) in the triangular plot defined in Fig. 4, for 
six representative values of Sm 2 /E u . The first subplot represents the averaged vacuum oscillation 
case. In the other subplots, the MSW effect is active and produces an increasing deformation of the 
curves as v e approaches ^3 (the left lower corner). For to = 0, tt/4, tt/2, however, the probability 
remains equal to the vacuum value, as discussed in the text. (The behavior at uj ~ cannot be 
graphically resolved in the last two subplots.) 
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3v MSW = P(y e -»v e ) 



10* 



10° 





10 2 




10 








1Z 


1 






-+-> 










10* 








io- 3 




10 L 




10 




10+ r 



— i i iiiiij i iiiiij ii imi| i ii ni| i 1 1 1 ■ i ■ | i ii i rrm^ — i i iiii^ i mi* 

6m a /E » 10"° eV'/MeV 




i iiiii| i i ■ ■ 1 1 ■ | i rrm^ — i i niiij i n ni| — i i iniij i i ini| i i mil] i i imij i i im 

. <5m ! /E= 1.75X10"° i 



Non-res. i Non-res. - 



0.9 






10 s 




10 










£Z 


1 


n 


-1 


-+-• 


10 








io- 2 




-3 




10 




-+ 
10 



1 10 10 1CT 10 10 10 10 10 10 10 1 10 10 1 
llj I Illllj llllllll I IIIIIJ I lllllll IIIIIJ I Illllj IIIIIJ I 



6mVE=1.62X10" 5 
Non-res. 



Non-res. -. r- 




I Illllj I 1 1 I Illllj lllllll| I Illllj lllllll| I Illllj lllllllj lllllll| I Mil 

, emVE=1.45X 10" 5 i 



Non-res. ! Non-res.-: 



0.9 



10 7 10 B 10 5 io -4 10 3 10 1 10 

1+ 






10 J 




10 5 




10 










c 


1 


o 






io" 1 








io" 2 




—3 




10 




io" 4 



iiiii| — i i iniij i iini| ii imi| i n ni| i 1 1 1 ■ i ■ | i ii iii^ i rrm^ — i i iiii^ i mi 

firn 2 /E=7.95XlO- 9 i 



0.9 



Non-res. j Non-res.' 

t 




I lllll| I I Illllj I lllll| — I I Illllj I II lll| — I I Illllj I I llllj I I lllll| I I Illllj I I MB 

. (5mVE:=5.15XlO- 7 




tci r\ Z M 



FIG. 7. As in Fig. 6, but in the bilogarithmic (tan 2 u, tan 2 (f>) representation. Notice the 
deformation of the iso-probability contours due to the A — ► □ mapping. In this representation 
the v±, V2, vj, mass eigenstates are reached respectively in the limits: tan 2 cj) — > and tan 2 cj — > 0; 
tan 2 and tan 2 tj — > oo; tan 2 <p — ► oo at any tan 2 tu. The thick, dashed lines represents the 

curves for maximal cj-mixing in vacuum and matter: to = it /A (vertical) and uj^ = ir/A (curve 
on the left). These two curves separate the zones where the neutrino propagation is non-resonant 
from those where it is resonant. 



32 



3v MSW = P(y e -»v e ) 



10 



I lllll| — I I lllllj I lllll| I I lllll| I II lllj I lllllj: 

tanV=0 




I iiiii^ I I lllllj I mnj — I I niiij I II uij — I I mpj i I i^iij 1 1 ijiiij j I I imij I 1 1 1 it 

- - - I = ! * 

Mm! 

D.6 




i niiij i iiiiiij i niiij i miiii| i iiiii| i inipi 
tanV=0.5 



> 



10 



10 



E 

<0 



10 



10 




I lllllj lllllllj I lllllj lllllllj I lllllj lllllllj I lllllj IJIIIIIj 

tanV = 1 I ! ! 



I lllllj I I 1 1 u 



0.4; 



Non-res. 




1 1 J 1 J 1 1 I 1 1 1 



J 1 J 1 1 tlll^ "' tillU i 1 1 

" 7 10" G 10" E 10"^ 10" S 10"^ 10" 1 1 10 10 J 10 3 10" 7 10" G 10" 5 10" + 10" S 10" 2 10"' 1 10 10 J 10 3 



— i i niiij i iiiiij ii imi| i ii wiy i iiii^i i n i iiiii| | i i iiii^ 

tanV= 1 ,5 




: 1 Mill] 1 Illlllj 1 lllllj 1 1 lllllj 


1 II lllj 1 1 lllllj 1 Mill 


l 


i imij 


1 1 II u 












Non-res. 




















: 






d> 






Res, 




CI 

O 






0.5 



















\an Z G) 



FIG. 8. Iso-lines of survival probability P(y e — > v e ) in the bilogarithmic (tan 2 lo, dm 2 / 1 E v ) plane, 
for representative values of tan 2 4>. The thick, dashed curves are defined as in Fig. 7. The resonant 
region has the characteristic triangular shape. The variations of P(y e — > u e ) in the plane decrease 
for increasing (p. 
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FIG. 9. Three-flavor MSW analysis of all solar neutrino data. The regions allowed at 90%, 
95%, 99% C.L. in the space {5m 2 , tan 2 cj, tan 2 <j>) are shown in planar (5m 2 , tan 2 uj) sections at 
twelve representative values of tan 2 4> ranging from to 1.5. The two- generation limit is recovered 
at tan 2 <j> = (first subplot). The two separate solutions at small 4> merge in one single solution 
for increasing (p. See the text for details. 
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FIG. 10. Values of A% 2 as a function of tan 2 <fi, for 5m 2 and tan 2 oj unconstrained. The value 
= 0, corresponding to the two-generation limit, is preferred, but values as high as tan 2 cp ~ 1.4 
(0 ~ 50°) are still allowed at 90% C.L. 
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FIG. 11. Three-flavor MSW analysis of individual neutrino data in the same coordinates as in 
Fig. 9, for three representative values of tan 2 cf>. The solid lines correspond to 95% C.L. contours. 
The allowed regions are marked by stars. Dotted lines represent iso-signal contours. 
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FIG. 12. Results of the three-flavor MSW analysis of all solar neutrino data. The regions allowed 
at 90%, 95%, 99% C.L. in the space (5m 2 , tan 2 a;, tan 2 <p) are shown in planar (tan 2 oj, tan 2 0) 
sections at twelve representative values of 5m 2 ranging from 1.5 x 10~ 4 to 2.2 x 1CT 6 eV 2 . The 
two-generation limit is recovered for tan 2 (lower side of each subplot). Notice the appearance 

of solutions disconnected from the two-generation limit. See the text for details. 
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FIG. 13. Results of the three-flavor MSW analysis of individual neutrino data in the same 
coordinates as in Fig. 12 for three representative values of 5m 2 . The solid lines correspond to 95% 
C.L. contours. The allowed regions are marked by stars. Dotted lines represent iso-signal contours. 
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